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A general  model  of  hlgh-energ^^  laser  interactions  with  solid  surfaces  is 
presented.  Fluid  transport  equations  are  used  to  describe  the  heating  and 
vaporization  of  a solid  surface  Irradiated  by  intense  laser  energy.  The  vapor- 
ized target  material  dlfhises  into  an  ambient  gas.  Both  the  target  vapor  and 
ambient  gas  can  ionize.  Separate  transport  equations  are  used  for  the  ambient 
gas,  the  gas  of  atoms  in  an  excited  state,  the  electrons,  the  target  vapor,  the 
ionlised  ambient  gas,  and  the  ionized  target  vapor.  Among  the  over  30  physical 
processes  Included  In  the  model  are;  laser  wave  absorption,  electron  and  Ion 
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20.  Abstract  /continued) 

gas  heating,  diffusion  of  each  species,  excitation  collisions,  recombination, 
radiation  transport,  photoionization,  shock  formation,  cascade  ionization, 
thermionic  emission,  neutral  impact  ionization  and  energy,  and  momentum 
transfer  among  all  fluid  species.  A computer  program  has  been  developed  to 
numerically  integrate  these  transport  equations.  The  ignition  and  propagation 
of  laser-supported  absorption  waves  (LAWs)  are  studied  as  a function  of  incident 
power  level.  Graphs  are  presented  of  the  temperature  and  density  profiles  of 
each  species  at  various  instants  in  time  as  a function  of  the  incident  laser 
power  level. 
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High-Energy  Laser — Target  Interactions 


1.  iNTRODUCTION 

The  interaction  of  high-energy  laser  radiation  with  a solid  surface  can  produce 
many  physical  processes.  These  processes  include  the  following:  absorption  of 
laser  radiation  by  a solid  target,  thermal  conduction  into  target,  vaporization  of 
target  material,  diffusion  of  target  material  into  an  ambient  gas,  thermionic  emis- 
sion from  heated  target  surface,  cascade  ionization  of  target  vapor  and  ambient 
gas  atoms,  heating  of  electrons  and  ions  by  inverse  bremsstrahlimg,  excitation 
and  de-excitation  collisions,  ener,  7 and  momentum  transfer  among  all  species, 
neutral  impact  ionization,  photoionization,  radiation  transport,  shock  formation, 
attenuation  of  laser  radiation  by  plasma,  diffusion  of  each  species,  electrical 
coupling  between  the  electrons  and  each  ion  species,  recombination,  electron 
impact  ionization  from  excited  states,  and  so  forth. 


(Received  for  publication  3 October  1975) 
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Previous  works  in  this  field  (Reedy,  ^ DeMtchella,  ^ Nielsen,  ^ Nielsen  end 

e e 

Genevan,  Harmon  ) have  Included  some  of  these  mechanisms,  but  not  ell  In  a 

self-consistent  theory.  Most  of  the  previous  studies  have  treated  the  problem 

using  a one-,  or  at  most,  two-fluid  model  with  separate  electron  and  heavy- 

7 R 

particle  teniperatures.  ' In  these  treatments,  the  ion  temperature  Is  assumed 
to  be  equal  to  the  heavy-particle  temperature.  However,  especially  In  the  Immed- 
iate vicinity  of  the  target  surface,  these  assumptions  may  not  be  entirely  valid, 
since  the  expanding  target  vapor  Is  much  hotter  than  the  ambient  gas.  Also,  since 
the  rate  of  diffusion  of  the  Ions  Is  coupled  to  that  of  the  electrons,  the  cooling  by 
expansion  for  the  Ions  will  be  different  from  that  for  the  neutral  species.  Finally, 
relaxation  times  between  species  can  greatly  Increase  so  that  the  time  for  energy 
equilibration  may  become  quite  large.  Indeed,  It  has  been  shown  that  by  assuming 
an  Initial  electron  density  of  about  10  part/cm  over  a distance  of  several  laser 

O 

wavelengths  (10.  6 pm),  f r Incident  laser  power  levels  greater  than  3 x 10  watts/ 
2 

cm  In  xenon,  the  Ion  temperature  can  be  different  than  the  heavy  particle  tem- 
perature. ^ 

P.  D.  Thomas has  performed  computer  solutions  for  the  ignition  of  laser 
absorption  waves  (LAW's)  baaed  on  a uniform  vaporization  model  that  takes  Into 
account  gas  dynamics  and  radiative  transport  in  the  interaction  region  over  the 
surface,  anu  for  heat  conduction  within  the  target.  Thomas^®  has  compr.red  his 


1.  Ready,  J.F.  (1971)  Effects  of  Hlgh-Power  Laser  Radiation,  Academic 

Press. 

2.  DeMlchells,  C.  (1969)  IEEE  J.  Quant.  Electronics  QE-5;188. 

3.  DeMlchells,  C.  (1970)  IEEE  J.  Quant.  Electronics  QE-6;630. 

4.  Nielsen,  P.  E.  (1974)  Breakdown  and  laser  absorption  waves.  Jour.  Defense 

Research.  Series  B,  December. 

5.  Nielsen,  P.  E.  and  Canavan,  G.  H.  (1974)  Theory  of  laser  target  Interaction, 

Jour.  Defense  Research,  Series  B,  December, 

6.  Harmon,  N.F,  (1974)  Proceedings  of  the  1973  DoD  Laser  Effects /Hardening 

Conference  (U),  Vol.  V,  Laser  supported  absorpllon  waves.  Mure  Corp.^ 
Rpt  M73-115. 

7.  Edwards,  A.,  Ferrlter,  N.M.,  Fleck,  J.A.,  and  Winslow,  A.M.  (1974) 

A theoretical  description  of  pulsed  laser-target  Interaction  in  an  air  environ- 
ment, published  in  Proceedings  of  the  1973  DoD  Laser  Effects /Hardening 
Conference  (U).  Vol.  V (■see  r ef.  6,  p.  115). 

8.  Stamm,  M.  R.  and  Nielsen,  P.  E.  (1974)  Nonequilibrium  effects  in  shock- and 

transport-induced  LAW  formation,  published  In  Proceedings  of  the  1973 
DoD  Laser  Effects /Hardening  Conference  (U).  Vol.  V (see  rer.  6,  p.  141). 

9.  Papa,  R.J.  (1974)  Propagation  of  Hlgh-Power  Laser  Radiation  In  Partially 

Ionized  Oases.  AFc3flL-TR-74-o632. 

10.  Thomas,  P.  D.  (1973)  Laser-Absorption  Wave  Initiation  In  Air  Over  Vaporiz- 
ing Targets  (s,ee  ref.  6,  p.  69). 
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calcuUtlotw  with  the  experlmentel  results  of  Msthemstlcsl  Sciences  Northwest^  ^ 

12 

end  Bsttelle  Memorlsl  Institute.  Themes  concluded  that  a uniform  vaporisation 

model  Is  probab?^  valid  at  10.  6*pm  radiation  for  Intensities  of  the  order  of  a few 
6 2 

times  10  watts /cm  or  leas,  but  the  physics  of  the  Interaction  Is  different  above 
this  Intensity. 

Musal^^  has  noted  that  in  the  laser  Intensity  range  3 x 10*^  to  3 x 10^  watts/ 

cm^.  laser-supported  absorption  waves  are  Ignited  In  times  much  shorter  than 

that  required  to  produce  uniform  surface  vaporisation.  Musal  has  presentsd  sn 

analytic  model  of  the  fleld-emlsalon/cascade-lonlaatlon  buildup  of  the  plasma  over 

a small  protrusion  on  the  laser-lrradlated  target.  Musal  has  concluded  that  the 

prompt  Initiation  of  laser-supported  absorption  waves  at  metallic  surfaces  at 

1 2 

Incident  laser-beam  Intensities  above  S x 10  watts/cm  a.:s  probably  caused  by 
cascade  air  breakdown  Initiated  from  localised  plasmas  generated  In  the  ambient 
atmosphere  by  Impact  Ionisation  from  electrons  field-emitted  from  surface  pro- 
trusions. 

Walters  has  Hated  the  following  mechanisms  for  LAW  Ignition: 

(1)  Uniform  target  vaporisation 

(2)  Oxide  vapor  absorption 

(3)  Chemical  reaction 

(4)  Thermionic  emission 

(5)  Planar  reflection-enhanced  fields 

(6)  Local  target  heating 

(7)  Contaminant  and  gas  desorption 
(d)  Field  emission 

(9)  Shock  heating  of  the  air 

(10)  Nonequilibrium  ionisation  of  the  vapor 

(11)  Defect  enhanced  fields 

Using  photographic,  photometric,  spectroscopic,  electronic,  and  surface  diagnos- 
tics, Walters  has  assessed  the  role  of  the  various  LAW  Ignition  mechanisms. 

He  has  concluded  that  thermionic  emission  and  local  target  heating  play  a major 
role  In  the  LAW  Ignition  process. 


11.  Klosterman,  E.  L. , Byron,  S.  R. , and  Newton,  J.F.  (1973)  Laser  Supported 

Combustion  Waves  Study  Report  No.  73-101-3,  Mathematical  Sciences 
itorthwesl,  btc.,  Seattle,  Wash. 

12.  Walters,  C.T.  (1972)  Review  of  Experiments  at  Battelle  presented  at  the 

Workshop  on  J^er-'Ehduced  Uomoustlon  Waves,  'Seattle,  Wash. , 

5-6  December  1972. 

13.  Muaal,  H.  M.  (1973)  Prompt  faltlatlon  of  Laser-Supported  Absorption  Waves 

In  Air  Via  Field-Emission  from  Target  Surfaces  (see  ref.  6.  p.  1571. 

14.  Walters,  C.T.  (1973)  Experlmwtal  Studies  of  LSD  Wave  Initiation  on 

Aluminum  Targets  (see  ref.  6,  p.  179). 
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In  tht  prMtnt  p«p«r,  m*chiinlam«  <1).  (3),  (4),  (6),  (B).  (8),  utd  (10>  «r« 
Included  in  e conetetent  feehlon.  Alto,  the  effect  of  the  Ionised  target  vapor  on 
the  electron  diffusion  process  is  properly  treated.  Cohen  et  al.^^  have  presented 
a method  of  analysis  by  which  the  electron  cascade  equations  are  coupled  to  the 
fluid  mechanics  conservation  equations  for  the  purpose  of  investigating  the  effect 
of  a shock-induced  flow  field  on  breakdown.  The  laaer  plasma  interaction  is 
treated  by  an  Idealised  single  fluid  model  which  includes  unequal  electron  and 
heavy  particle  temperatures  and  a two-stage  ionisation  model.  A similar  model 
has  been  used  by  Stamm  and  Nielsen.^  However,  the  transport  equations  preoentsd 
1.1  references  8,  10,  and  15  are  such  that  some  of  the  nonlinear  convective  deriva- 
tive terms  are  neglected.  Also,  in  the  equations  for  the  time  rate  of  change  of 
energy  density,  some  of  the  terms  representiog  the  time  rate  of  change  of  the 
fluid  energy  due  to  its  convective  motion  are  neglected.  At  high  laser  intensities, 
these  terms  may  not  be  negligible. 

Part  of  the  purpose  of  the  present  report  is  to  relax  some  of  the  assumptions 
that  were  made  in  references  1,  7,  8,  10,  13.  and  15,  in  order  to  Investigate  their 
range  of  validity.  For  example,  it  may  not  always  be  appropriate  to  write  down 
the  continuity  and  energy*  equations  in  a single  LAgrangian  frame  for  all  species, 
because  this  Implies  that  the  coi.vective  velocities  of  electrons,  ions,  and  neutrala 
are  all  equal.  Also,  the  electron  diffusion  coefficients,  which  are  normally  de- 
rived in  an  Eulerian  (fixed  frame),  will  differ  in  a LAgrangian  frame,  m addition, 
since  the  laser  flux  levels  can  reach  quite  high  values,  the  nonlinear  convective 
derivative  terms  in  the  n'omentum  transfer  equation  for  the  electrons  should  not 
be  neglected  a priori.  However,  it  can  be  shown  that  for  10.  S-|im  flux  levels  up 
to  5 X 10  watts /cm  , the  nonlinear  convective  derivative  terms  in  the  electron 
gas  can  Indeed  be  neglected  (see  Papa,^  p.  54).  Finally  since  the  target  ion  density 
may  become  quite  large  in  the  vicinity  of  the  target,  the  effect  of  the  presence  of 
two  distinct  iwi  species  (target  and  ambient  gas)  on  the  electron  diffusion  coeffi- 
cient should  be  investigated.  Ine  fluid  transport  equations,  written  in  an  Eulerian 

1C 

frame  with  no  terms  omitted,  are  given  in  Holt  and  Haskell,  pp.  156  to  166. 

These  general  equations  are  the  ones  used  in  this  report. 

The  present  report  Includes  the  following  physical  mechanisms  describing 
high-energy  laser -target  interactions: 

(1)  lAser  wave  absorption  in  a plasma  and  in  the  target; 


15.  Cohen,  H.D.,  Su,  F.  Y.,  and  Boni,  A. A.  (19*^3)  A Simple  Nonequiltbrlum 

Model  for  a Shock-Heat^  Monatomic  Gas  IrradtaiedW  an  mtwse  LAser 
Beam  (see  ref.  6.  o.  ^^9). 

16.  Holt.  E.H.  and  .'laskell,  R.E.  (1965)  Foundations  of  Plasma  Dynamics, 

Macmillan. 
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(2)  Electron  MVd  ton  gn«  heating  by  tnveree  brenteetrehlung  due  to  absorption 

of  laser  wave  energyt 

(3)  Electron  diffusion  coupled  with  the  dlffbalon  of  the  ambient  gaa  Iona  and 
the  target  vapor  lonet 

(4)  Excitation  of  one  energy  level  by  electron  Impact  on  a neutral  atom: 

<S)  De-excttatlon  of  energy  level  by  radiation: 

(6)  Diffusion  of  resonant  radiation  and  subsequent  excitation  of  neutrals  by 
resonant  radiation  absorption: 

(t)  Electron  Impact  ionisation  from  ground  state: 

(21  Electron  impact  lonlaatlon  from  excited  level: 

(9)  Diffusion  of  neutral-state  atoms  and  excited-state  atoms  dut)  to  tempera- 
ture and  density  gradients: 

(101  Recombination: 

(a)  radiative, 

(b)  electron-electron-l<m, 

(c)  transitional: 

(11)  Momentum  transfer  from  laser  wave  to  electron  gaa  due  to  Inhomo- 
genltlea  in  dielectric  coefficient  (nonlinear  forces  discussed  by  Hora^*^): 

(12)  Thermal  conduction  in  electron  gas,  ion  gas,  neutral  gaa,  gas  of  excited 
atonw  and  target  vapor: 

(13)  Bremsstrahlung  loss  by  electrons  colliding  with  ions: 

(14)  Radiative  recombination  cooling  of  electron  gas: 

(15)  Energy  transfer  by  collisions  between  electrons  and  ions: 

(12)  Energy  transfer  by  colllstona  between  electrons  and  neutrals: 

(17)  Energy  transfer  by  collisions  between  Ions  and  neutrals: 

(16)  Nonlinear  convective  derivative  terms  are  retained  in  the  momentum 
transfer  equations. 

(19)  Heating  of  the  ions  due  to  dynamic  polarisation  of  the  plasma  (Vinogradov 
and  Pustovalov**): 

(20)  Arbitrary  laser  frequency  (so  long  as  oc  where 

^excited  ‘ •***'^6y  excited  state,  h * Planck's  constant  and  * laser  frequency, 
so  that  classical  physics  applies): 

(21)  Arbltrsry  pulse  length,  pulse  shape,  and  pulse  repetition  frequency: 

(22)  Shock  formation  mechanisms: 

(23)  Photo-lonlsatlon  due  to  continuum  radiation  absorption  from  the  ground 
state  or  from  an  excited  state: 


17.  Hora,  H.  (1989)  Fhya.  Fluids  12:182. 

18.  Vinogradov,  A.V.  and  Pustovalov,  V.  V.  (1972)  Soviet  Jour.  Quant. 

Electronics  2:91. 
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(24)  Neutral-neutral  impact  ionization; 

(25)  Absorption  of  laser  energy  into  a target  and  thermal  conduction  into  the 
target; 

(26)  Vaporization  of  target  material,  with  diffusion  and  ionization  of  target 
material  into  an  ambient  gas; 

(27)  Thermionic  emission  from  target  surface; 

(28)  Effect  of  target  ions  on  electron  ambipolar  diffusion  coefficient; 

(29)  Energy  and  momentum  exchange  between  target  vapor  and  ambient  gas 
particles; 

(30)  Energy  and  momentum  exchange  between  target  vapor  and  electrons; 

(31)  Ionization  of  target  vapor  by  electron  impact; 

(32)  Ionization  of  target  vapor  by  neutral -neutral  Impact. 

9 

In  a previous  report,  it  was  explained  how  multi-fluid  transport  equations 
can  be  used  to  describe  the  physical  mechanisms  (1)  through  (24).  Mechanism  (25), 
absorption  of  laser  energy  into  a target,  and  thermal  conduction  into  the  target, 
are  described  in  this  present  report  by  solving  the  one -dimensional  heat-flow 
equation.  The  absorptivity  and  reflectivity  of  the  target  are  functions  of  the  inci- 
dent laser-power  level.  In  Fig,  3.26,  p.  116  of  reference  1,  the  reflectivity  of 
various  materials  is  plotted  as  a function  of  laser-power  level  for  a Q-switched 
Nd-glass  laser.  In  the  present  report,  it  is  assumed  that  the  reflectivity  of  a 
target  has  the  same  behavior  under  Irradiation  at  10.  6 pm  wavelength.  The 
solution  to  the  heat -conduction  equation  for  an  arbitrary  laser-pulse  shape  is 
given  by  Eq.  (3.9),  p.  73  of  reference  1: 


1 /2  /• 

T,j,(x,  t)  = i J dr  exp  (-x^/(4Kt)) 

where 

T,p  = temperature  of  target; 

X = spatial  position,  with  x = 0 taken  as  target  surface; 
t = elapsed  time; 

K = thermal  dlffU8l”  ty; 

K = thermal  condv.  .tlvity; 

F = absorbed  laser  flux  as  a function  of  time. 


(1) 
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The  target  ia  aasumed  to  vaporize  according  to  a model  proposed  by 
19 

Anisimov.  The  surface  temperature  of  the  target  is  obtained  by  equating  the 
absorbed  laser  energy  density  to  the  surface  temperature  rise  times  the  specific 
heat  plus  the  specific  heat  of  vaporization: 

F = i'ssPtS’^' , (2) 

where 

V = velocity  of  retreating  target  surface, 
ss 

p = target  density, 

Cp  = heat  capacity  per  unit  mass, 

T'  = maximum  target  surface  temperature, 

L = specific  heat  of  vaporization  per  unit  mass. 

If  N,p(z  = L)  is  the  density  of  target  vapor  atoms  at  the  target  surface,  then 
dN„(L) 

‘'ss  = 


where  h^  = thickness  of  an  atomic  layer  In  target  material. 

The  continuity  equation  for  the  target  vapor  just  at  the  target  surface  is  given  by 


exp  (-LM/N^kT') 


where 

= Debye  frequency, 

M = atomic  weight  of  target  material, 

= Avogadro's  number, 
k = Boltzmann's  constant. 

Combining  equations  (2),  (3),  and  (4)  gives  an  implicit  equation  for  the  maximum 


target  surface  temperature  T': 

F = ha  VoP  [CpT'  + L]  exp  (-LM/N^kT')  . 


19.  Anisimov,  b.  i.  '1967)  Soviet  Phys.  Tech.  Phys.  11:945. 
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The  fluid  transport  equations  governing  the  target  vapor  density  and  tempera- 
ture are  similar  to  the  transport  equations  for  the  other  species,  such  as  ambient- 
gas  atoms,  excited-state  atoms,  electrons,  ambient-gas  Ions,  and  target  vapor 
Ions  (Holt  and  Haskell,  pp.  156-166). 

A very  general  computer  program  has  been  developed  to  solve  this  complica- 
ted set  of  nonlinear,  partial  differential  equations.  The  program  gives  the  laser 
Intensity  at  any  point  from  x = -«  up  to  the  target  surface  x = L.  The  attenuation 
of  the  laser  beam  due  to  any  Ionized  particles  generated  by  the  interactions  Is 
accounted  for  by  representing  the  laser  electric -field  Intensity  through  the  WKB 
approximation.  Equation  (1)  is  solved  numerically  for  the  target  surface  tempera- 
ture, Tij(x  = L).  Then  Eq.  (2)  la  solved  for  T'  using  the  Newton-Raphson  method 
(Stark, p.  85);  but  T,p(x  = L)  cannot  exceed  T'.  Equation  (4)  is  solved  for  N,j„ 
the  target-vapor  density,  just  at  the  target  surface.  The  fluid  transport  equations 
for  the  target-vapor  density  and  temperature,  ambient-gas  density  and  tempera- 
ture, excited -state  atom  density,  electron  density  and  temperature,  ambient -gas 
ion  density,  target-vapor  ion  density,  and  ion  temperature  are  coupled  to  Maxwell's 
equations  which  govern  the  laser  wave-propagation  characteristics. 

As  explained  in  a previous  report,  the  plasma  transport  equations  are  of  the 
telegraph  type  (Holt  and  Haskell,^®  p.  186).  The  coupled,  partial  differential 
equations  are  converted  into  a system  of  difference  equations.  A partially  implicit, 
partially  explicit  (PI-PE)  scheme  for  numerically  integrating  these  transport  equa- 
tions is  used  for  purposes  of  maintaining  numerical  stability,  even  when  a large 
time  step  is  selected.  The  PI-PE  technique  for  integrating  transport  equations  in 
time  is  discussed  by  Rlchtmyer  and  Morton.  The  resulting  system  of  difference 
equations  for  the  particle  densities,  flaxes,  and  temperatures  is  solved  at  each 
time  step  by  inverting  a band-structured  matrix.  An  efficient  computer  algorithm 
has  been  developed  for  inverting  such  band-structured  matrices. 

At  each  time  step,  the  electromagnetic-field  distribution  and  the  transmission 
coefficient  are  calculated  by  using  the  WKB  method  to  solve  the  wave  equation. 

Also,  if  the  percent  change  of  all  species  is  less  than  a certain  tolerance  from  one 
time  step  to  the  next,  the  time  step  is  increased.  If  the  percent  change  of  any  one 
species  is  greater  than  a certain  tolerance  from  one  time  step  to  the  next,  the  time 
step  is  decreased.  Finally,  if  the  percent  change  of  density  or  temperature  of  any 
one  species  is  less  than  a certain  tolerance  for  ten  consecutive  time  steps,  that 
particular  species  is  "frozen  out"  for  ninety  consecutive  time  steps. 

The  various  species  that  are  involved  in  the  high-energy  laser  — target  inter- 
actions include:  electrons,  target-vapor  ions,  ambient-gas  ions,  ambient  gas  atoms 

20.  Stark,  P.  A.  (1970)  Introduction  to  Numerical  Methods.  Macmillan. 

21.  Rlc.htmyer,  R.D.  and  Morton,  K.W.  (1967)  Difference  Methods  for  Initial 

Value  Problems,  Interscience.  
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excited -state  atoms,  and  target  vapor.  Thus,  there  are  six  sets  of  coupled  fluid 
equations.  For  a given  set  of. Initial  conditions,  laser-pulse  shape,  pulse  height, 
pulse  length  and  target  conditions,  graphs  are  presented  of  the  particle  densities 
and  temperatures  as  a function  of  position  at  different  Instants  in  time.  In  this 
particular  report  the  ambient  gas  is  chosen  to  be  xenon,  the  laser  wavelength  la 
10.  6 /sm,  and  the  target  Is  aluminum. 


2.  MULTIFLUID  TRANSPORT  EQUATIONS 

In  reference  9,  the  multlfluld  transport  equations  were  given  for  a high- 
Intenslty  laser  wave  propagating  through  a partially  Ionized  plasma.  These  equa- 
tions were  presented  after  making  27  self-consistent  assumptions.  In  the  present 
report,  target  heating  due  to  laser  absorption,  target  vaporization,  diffusion  of 
target  vapor  and  Ionization  of  target  vapor  are  described,  so  that  Ignition  and 
propagation  of  LAW'S  may  be  systematically  Investigated.  Thus,  in  addition  to 

Q 

the  assumptions  listed  In  the  report  by  Papa,  several  additional  assumptions  are 
added  In  order  to  Include  the  target  Interactions. 

A list  of  all  these  assumptions  follows: 

(1)  It  Is  assumed  that  there  are  no  variations  of  the  plasma  parameters  in  the 
y or  z direction,  only  in  the  x direction  (direction  of  laser  wave  propagation).  The 
variation  of  the  electric  field  intensity  Is  also  neglected  in  the  y and  z directions. 
These  assumptions  are  reasonable,  provided  the  diameter  of  the  laser  beam  is 
much  larger  than  any  characteristic  lengths  in  the  plasma  (such  as  Debye  length, 
diffusion  length,  and  scale  lengths  for  temperature  and  density  gradients). 

(2)  The  heat  flux  tensor  is  assumed  to  be  Isotropic  for  thermal  conduction  In 
the  electron  gas,  ion  gas,  and  neutral  gas. 

(3)  There  is  only  one  excitation  potential  of  the  neutral  gas. 

(4)  The  electric  field  vector  E^.  is  transverse  to  the  direction  of  wave  propa- 
gatioi..  The  incident  laser  wave  is  a plane  wave. 

(5)  There  are  no  inelastic  collisions  of  the  second  kind;  only  electron  impact 
and  absorption  of  resonant  radiation  can  excite  the  neutral  particles  to  the  energy 
level  Egjj. 

(6)  There  is  no  multiphoton  absorption. 

(7)  The  electron  and  ion  energy  changes  are  assumed  to  occur  continuously, 
so  that  the  electron  energy  gain  and  loss  terms  are  treated  according  to  the  laws 
of  classical  physics.  This  assumes  the  laser  photon  energy  is  1 eV  or  lejs  (CO» 

A 

laser). 

(8)  It  is  assumed  that  at  t = 0 the  electron-density  profile  N^(x,.t  = 0),  the 

ion -density  profile  Nj(x,  t = 0),  and  the  neutral-gas -density  profile  Njj(x,  t = 0)  are 
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known.  Also,  ths  electron  temperature  T^(x,  t = 0),  ion  temperature  Tj(x,  t = 0) 

and  neutral-gas  temperature  are  assumed  to  be  known  at  t = 0.  These  initial 

electrons  and  ions  are  not  due  to  the  prebreakdown  ionization  observed  by 
22  

Papoular,  since  LAW'S  may  be  ignited  from  targets  at  laser  fluxes  lower  than 
that  required  to  initiate  breakdown  in  clean  gases  without  a target.  These  initial 
conditions  may  have  been  the  result  of  several  different  types  of  phenomena,  such 
as  the  direct  emission  of  particles  from  the  target  surface  due  to  either  photo- 
electric effect,  thermionic  emission,  multiphoton  absorption,  or  the  tunneling 
23 

effect.  Keldysh  has  shown  that  a frequency-dependent  tunneling  mechanism 
becomes  equivalent  to  multiphoton  ionization  in  the  high-frequency  limit. 

(9)  The  possibility  of  shock  generation  occurring  in  the  emerging  target  vapor 
will  be  neglected:  shocks  can  form  only  in  the  ambient  gas. 

(10)  Each  atom  is  assumed  to  be  capable  of  being  only  singly  ionized  so  that 
tnultiple  ionization  effects  are  neglected.  Only  the  following  ionization  mechanisms 
are  considered: 

(a)  neutral -neutral  gas  atom  impact  ionization 
Xe  + Xe  — Xe^  + Xe  + e ; 

(b)  neutral-neutral  target  vapor  atom  impact  ionization 
A1  + A1  -•  A1  + Al”*"  + e : 

(c)  neutral-target  vapor  atom— neutral -gas  atom  impact  ionization 

A1  + Xe  — Xe'^  + A1  + e ; - 

(d)  electron-neutral  gas  atom  impact  ionization 
e + Xe  -•  Xe^  + 2e  ; 

(e)  electron-neutral  target  vapor  atom  impact  ionization 
e + A1  -►  Al"*^  + 2e  ; 


22.  Papoular,  R.  (1972)  The  initial  stages  of  laser  induced  gas  breakdown, 

article  in  Laser  Interaction  and  Related  Plasma  Phenomena,  Vol.  2, 
Proc,  of  the  Second  Workshop  held  at  RPl,  Hartford  Graduate  Center. 
edited  by  Schwarz  and  H Hora,  Plenum  t’ress. 

23.  Keldysh,  L.  V.  (1965)  Soviet  Physics  JETP  21:1135. 
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(f)  photo-lonizatlon  of  the  neutral-gas  atoms  by  radiation  from  free- 
free  transitions 

hv  + Xe  -►  Xe^  + e s 

(g)  photolonlzatlon  of  excited-state  atoms  by  radiation  from  free-free 
transitions 

hw  + Xe'  — Xe'*’  + e . 

(11)  Inelastic  Ion -neutral  collisions  are  neglected,  which  Implies  that  charge- 
transfer  reactions  and  Ion-neutral  Impact  Ionizations  are  neglected. 

(12)  The  Thomson  cross  section  for  scattering  of  laser  radiation  by  free 
electrons  Is  assumed  to  be  negligible. 

(13)  Light  scattering  from  electron-density  fluctuations  as  discussed  by 
24 

Ramsden  la  negligible. 

(14)  The  heating  cf  the  electrons  by  stimulated  Compton  scattering,  as  dls- 

18 

cussed  by  Vinogradov  and  Pustovalov,  la  negligible.  However,  the  heating  of 
the  Ion  gas,  due  to  the  dynamic  polarization  of  the  plasma  at  the  flrequc-ncy  of  the 
scattered  photons,  Is  Included. 

(15)  The  radiation  at  frequency  = E^jj/h  la  assumed  to  have  a very  small 

mean-free  path  In  the  plasma.  Thus,  the  plasma  is  optically  thick  at  this  frequency 

and  only  the  wings  of  the  resonant  radiation  line  can  diffuse.  ThU  Implies  that  the 

boimd-bound  absorption  cross  section  is  large.  The  diffusion  of  the  resonant  radl- 

25 

ation  Is  treated  In  the  manner  discussed  by  Myshenkov  and  Raider. 

(16)  The  electron-dlstrlbutlcn  function,  the  lon-dlstrlbutlon  function,  the 
neutral  ambient-gas  atom  distribution  function,  and  the  target -vapor  atom  dis- 
tribution function  are  all  assumed  to  be  close  to  Maxwellian,  so  that  one  may 
define  an  effective  electron  temperature  T.,  an  effective  Ion  temperature  T.,  an 
effective  neutral  gas  temperature  and  a target  x'spor  temperature  T.^.  Thus, 
the  Maxwelllanlzatlon  times  for  each  species  are  assumed  to  be  shorter  than  any 
time  scales  for  temperature  or  density  changes. 

(17)  The  neutral  gas  particles  In  the  excited  state  are  assumed  to  have 
the  same  translational  temperature  T^  as  the  atoms  In  the  ground  state. 

(18)  The  electron  pressure  tensor,  ths  Ion  pressure  tensor,  the  neutral  gas 
pressure  tensor,  and  the  target  vapor-pressure  tensor  are  all  assumed  to  be  Iso- 
tropic. 

24.  Ramsden,  S.A.  (1968)  Physics  of  Hot  Plasmas.  Scottish  Universities  Summer 

School,  Plenum  Press. 

25.  Myshenkov,  V.  I.  and  Ralzer,  Yu.  P.  (1972)  Soviet  Phys.  JETP  34:1001, 
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(19)  There  Is  no  dc  magnetic  field  present. 

(20)  There  Is  no  effect  on  the  plasma  processes  due  to  metastable  states. 

(21)  The  plasma  Is  assumed  to  remain  nonrelatlvlstlc,  so  that  the  effect  of 
the  ac  magnetic  field  Is  negligible  compared  wU*'  the  ac  electric  field  (except 
Insofar  as  the  ac  magnetic  field  contributes  to  the  nonlinear  force  discussed  by 
Hora^^). 

(22)  The  wave  period  T = 2r/u  of  the  laser  radiation  Is  much  shorter  than  the 

time  scale  for  electron  density  and  temperature  changes,  so  that  the  electric  field 

26  27 

Ey  remains  monochromatic  In  the  time-varying  plasma  (Papa  * ). 

Ey  oc  exp  (-iut). 

(23)  Electron  attachment  1s  neglected  and  It  Is  assumed  that  there  are  no 
negative  ions  present  In  the  plasma. 

(24)  The  thermal  DeBroglle  wavelength  of  the  electrons  much 

less  than  the  average  particle  spacing 


DeBroglle 


« N 


-1/3 


e 


where 

Ng  = electron  particle  density, 

^DeBroglle  ' *‘/(2irmkTg)^/^  , 
h = Planck's  constant, 
m 3 electron  mass, 
k = Boltzmann's  constant, 

Tg  = electron  temperature. 

Then,  the  laws  of  classical  physics  can  be  used  to  describe  the  plasma. 

(25)  The  nonlinear  force  term 


(6) 


17 

discussed  by  Hora  affects  only  the  electron  momentum  equation,  but  not  the  Ion 
equation, 

(26)  The  term  fj^j^  • v^^  Is  neglected  In  the  energy  equation  for  the  electrons 
compared  with  the  joule  heating  term  • Ey,  where  Jy  = ac  current  density,  and 
Vjj  = electron  velocity  In  the  x direction. 


26.  Papa,  R.J.  (1965)  Can.  J.  Physics  43;38. 

27.  Papa,  R.J.  (1968)  Can.  J.  Physics.  46:889. 


(27)  Th«  target  is  assumed  not  to  melt,  the  target  material  passes  directly 
from  the  solid  to  the  vapor  phase. 

(28)  Photo 'ionization  occurs  only  in  the  ambient  gas,  not  in  the  target  vapor. 

(29)  The  nonlinear  convective  derivative  terms  are  neglected  in  the  momen- 
tum transfer  equations  for  the  elections,  target  vapor  ions,  ambient  gas  ions  and 
target  vapor  atoms . This  permits  the  particle  fluxes  to  be  written  in  terms  of 
appropriate  difhislon  coefficients  multiplied  by  density  and  temperature  gradients. 

(30)  The  laser  wave  that  is  reflected  from  the  target  surface  is  neglected, 
there  is  only  a wave  travelling  toward  the  target. 

(31)  The  ambient  gas  ions  are  assumed  to  have  the  same  temperature  T^  as 
the  ionized  target  vapor. 

(32)  The  ac  current  density  is  made  up  of  a component  due  to  electron 

A ^ I 

motion  and  to  the  motion  of  the  ambient  gas  ions  ; the  current  density  due 

to  the  motion  of  the  heavy  target  vapor  ions  is  negligible: 


With  these  assumptions,  Maxwell's  equations  coupled  with  the  fluid  transport 

equations  for  the  electrons,  ions,  neutral  gas  atoms,  excited  state  atoms,  target 

1 6 

vapor  ions,  and  neutral  target  vapor  become  (Holt  and  Haskell,  p.  156  et  seq. ): 
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+ (qNj,  + q„N^,)(3kT,/2) 

■*■  [‘'nITI^N^*^!  '‘'  '*’  ‘'nNIT^t]  <3kT^/2)  . 


(V) 


*^ex‘'exl^^l/’^"’l^  ■ ‘'exl^ex  ‘‘‘  \x‘'exN  ' ‘'exN^e 

■•■  ^ex*'exNT^T/^T^  ■ ‘'exNT^ex  * 


(19) 


(20) 


t 


1 

I 

i 


[ 

1 

) 

) 
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■^'iTi'efi)W~  il~  nT~  ^1^1 1#  iti  ~ ^'lyhVTi  ft  T~f  i ''  Ji~t  i^i^rirr  wifai-''i*~'iifiirt  ii  Viiiii«*vi  frcTfiin  f~ai 


rfjgja-fcvM i V .: 


It  (t"!",)  *1  »|i  iT,n'  *T”ilt  *T"'ili<'iK> 


8T 


-fe«i  (»,  ■ Oav  ‘ “.iNil" ■ V '.i 


■ •’'i  - V ■'m  -r«.  ••'i-'  "it’  {4‘’’i) 

*'nnit**t]  <aaTj/a> 

* (it)  (-^)  (;;V)  (ifc)  •'  ‘l*’  (-^) 

X N (J . a.a . s.*) 

® (T^  + Tjr  [ <2  + a*)  V^/ 


ti/a 


X tan“^  (^/«) 


(21) 


o ®(NrpT»p)  A w * F|» 

l^dt" 


■ '*’t^  *'eNT  ^i'*'  °NNT*^N^'^N  ' ‘'nnT 

- (lr*‘'^T)'^®r^e”rr  (t'^'^t) 

■ft  **  lx  ■*■  I *"NTt  * 7 “nTx 


(22) 


8T 


Tx-N^=  QeNNel  <T,  - V ^ VA  <^i  ' V v 


IN 


'*'  ®NNT  **N  *'nNT  ' ''l**e**N 
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- ■'sm''N<'*T  * V*>  - "h  (i^^N) 

♦ N, N,  Mq  Nj,  + N„J  (S  kT„/a)  . (23  contd) 


8B 


®»E  * I 


8(N.  V •) 


(24) 

(25) 

(26) 
(27) 


»(Niv3 

gt~^ 


(N,t,‘v^‘)  - N, 


-"i-mV 


-N,»U<V 


-v') 

y 


(28) 
. (28) 


Her«,  the  leaer  wave  propagates 'm  the  + x direction  and  the  target  la  located  at 
X - L. 


E„  and  H oc  exp  (-Iw  t)  . 

y * 

Nj  < particle  number  denalty. 

j * e -*  electrons. 

> I -•  ambient  gas  Iona. 

- IT  — target  vapor  Iona. 

> N -*  ambient  gas  atoms. 

X NT  — target  vapor  atoms. 

> ex  ->  excited  state  atoms. 

> nux  density  of  species  j. 
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-r.-  .. 


“^NNIT 


■ ttmperfttur*  of  apocloa  J. 

* lonlBatlon  froqutncy  du*  to  oloctron  Impact  on  ambient  gaa  atom# 
(iMitaation  croea  aaetton,  timaa  ralatlva  valoetty). 

> tontiatton  fraquancy  dua  to  alactron  Impact  on  target  vapor  atoma. 

« loniaatlon  fraquancy  due  to  a neutral  atom  colliding  with  an  ambient 
gaa  atom. 

« loniaatlon  frequency  dua  to  a neutral  atom  colliding  with  a target 
vapoi  atom. 

> loniaatlon  fraquancy  dua  to  alactron  Impact  on  an  excited  atata  atom. 

> total  recombination  rata  coefficient  for  radiative  recombination  and 
alactron-alactron-lon  recombination. 

* photo-lonlaatlon  rata  for  loniaatlon  of  ambient  gaa  atoma  by  continuum 
radiation. 

■ photo-lonlaatlon  rata  for  loniaatlon  of  excited  state  atoma  by  con- 
tinuum radiation. 

* diffusion  coefficient,  for  resonant  radiation  (bound-bound  tranaltlona) 

« (I*/3t). 

■ mean  free  path  of  reaonant  radiation  quanta. 

> average  lifetime  of  excited  atomic  atate. 

■ rate  coefficient  for  production  of  excited  states  by  electron  Impact 
on  ambient  gaa  atoms. 

= average  time  for  resonant  radiation  to  escape  from  plasma  (laterally, 
In  y-s  plane). 

= Boltsmann's  constant. 

X particle  mass  of  species 

■ electron  mass, 

« ac  electric  field  of  laser  wave. 

a ac  magnetic  field  of  laser  wave, 

> particle  density  of  target  vapor. 

« charge  on  electron. 

> particle  velocity  in  y direction  of  species  j. 

= particle  velocity  In  x direction  of  species  j. 


•o 

Mo 

''J 

** 

"tN 

‘'.I 


c 


h 

E 

E 

E 


1 

IT 

ex 


Qt 


■ permUttvtty  of  free  space. 

• parmeabllUy  of  free  space. 


■ quasl-statle  electric  field  due  to  charge  aeparatlon. 

* electron*neutral  collision  frequency  for  momentum  transfer. 

> electron-ion  collision  frequMcy. 

■ eoUlalon  frequency  for  momentum  transfer  between  species  k and 
species  j. 

* nij). 

■ collision  frequency  for  momentum  transfer  between  species  k end 
species  ] transformed  from  center  of  mass  coordinates  to  laboratory 
coordinates. 

- thermal  conductivity  of  gaa  consisting  of  species  J. 

« fractional  energy  loss  In  collisions  between  speclea  k and  species 
2 

j ■ 2m|^mj/(m|^  + m^)  for  elastic  collisions  between  species  k and 
species  ]. 

* velocity  of  light. 

■ Planck's  constant. 

= Ionisation  energy  of  ambient  gas  from  ground  state. 

= Ionisation  energy  of  target  material. 

> energy  of  excited  atomic  state. 

> rate  of  emission  of  target  vapor  at  target  surface, 

* exp  (-LM/Nj,kT^(x  = L)1 . 

= Debye  fiequency. 

* atomic  weight  of  target  material. 

= Avogadro's  number. 

: diffusion  coefficient  for  target  vapor, 

* ac  current  density  of  species  k. 

* e|Ey|/(mu). 

= frequency  of  laser  radlatlcm. 
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L 

2 

P 


o 

2 


2" 

= - u_  . 

p 

2 

= square  of  plasma  frequency  = N e /(me_). 

o o 

= relative  frequency  spread  of  laser  radiation. 
= (me^/8k)(Au/u)^. 

= i(mc2/kTg)(u2/<u,2 


The  physical  significance  of  each  term  in  Eqs.  (7)  through  (29)  is  discussed  in 
reference  (9).  In  reference  9,  it  was  explained  how  additional  assumptions  are 
made  to  further  simplify  Eqs.  (7)  through  (29).  The  terms 

V e §. 

'^y  3t 


and 


a 

■55T 


in  Eq.  (28)  for  the  y component  of  the  electron  velocity  are  small  and  can  be 
neglected  if  it  is  assumed  that; 

(33)  The  time  scale  for  electron-density  changes  is  long  compared  with 
a wave  period  ® 


Tj^  55>  (27r/u) 


(34)  The  electron-neutral  collision  frequency  or  the  electron- ion  collision 
frequency  is  large  compared  with  , where  is  the  space  scale  for 

electron  density  gradients  ® ® 


"eN^  /L 


N, 


or 


■'.I*  V 'l- 


N. 


Similar  assumptions  can  be  made  in  Eq.  (29)  for  the  y component  of  the  ion  velocity: 

(35)  It  is  assumed  that  the  time  scale  for  ion-density  changes  Is  long  com- 
pared with  a wave  ; -iriod  ^ 

Tj^  •»  (2x/u)  . 


'I 
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(36)  It  Is  assumed  that  the  Ion -neutral  collision  frequency  Vjjj  or  the  Ion- 
electron  collision  frequency  be  large  compared  with  where  Is  the 

space  scale  for  Ion-density  changes 


or 


With  assumptions  (33)  through  (35),  Eqs.  (28)  and  (29)  become: 

■ ■‘'cnV  * ■'el  - V'’  • 

-•“V  - (Sf)  -y 

Solving  Eqs.  (30)  and  (31)  yields: 


(30) 

(31) 


e +eE^<‘“-7j^) 

'^y  - m : DENCKT 


(32) 


I ■ *'eN^ 

'^y  " m^  » DETJOir 


(33) 


where 


DENOM  = -0)2  . lu  (v^j,  + + vj^)  + ‘'eN^Ie  ‘'el  ‘'iN  ' 

The  self-consistent  electric  field  that  appears  In  Eqs.  (13),  (14),  (15),  and 
(19)  may  be  eliminated  by  making  several  additional  assumptions.  The  Initial 
Ionization  is  assumed  to  be  sufficiently  high  that  the  Debye  length  f p la  much  less 
than  the  scale  lengths  and  fo^^ectron  and  ion  density  changes: 

(37) 


where 


-2 


'd 


/2,  nFn  <N,*N,t|1 

■W°)K  J ■ 
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From  Poisson's  £q.  (19)  and  assumption  (37)  it  is  easily  shown  that 

Ne  = Nj  + . C 

for  distances  greater  than  a Debye  length.  Then,  subti  acting  Eqs.  (9)  and  (10) 
from  Eq.  (7)  and  using  Eq.  (34)  gives: 

8x  9x  8x  ■ '■ 

From  assumption  (19)  one  may  deduce  that  the  plasma  is  irrotational,  so  that 
Eq.  (35)  implies  that 


r = r + r 
e ^ I ^ IT  ‘ 


In  order  to  simplify  Eqs.  (13),  (14),  and  (15),  the  self-consistent  electric 
field  E^  may  be  eliminated  by  making  the  following  assumptions: 


(38) 

9r 

— ^<<.  V r 
at  eN  e 

(39) 

T « • 

el  eN 

(40) 

*^eex  ^ *^eN  ' 

(41) 

‘'eNT  ‘'eN  ' 

(42) 

(43) 

7,  « 7,.,  . 
le  IN 

(44) 

7t  « 7^._  . 

lex  IN 

(45) 

"iNT  "'in  ■ 

(46) 

OF  IT 

at  « "ITN^ 

(47) 

"iTe  "iTN  ■ 

(48) 

"lex  "iTN  • 

r_  = -DIFFE,  - DIFFE, 


8^ 


-ANLP|j[-i-tjEl2+|^ol«l^]  • 


= - DIFFITj  - DIFFITj  ^ - ANLPI  j^-g-e^l  e|  ^ + y^i^l  h|  ; 


where 


DIFFE 


.■  (Dk)(-^)[^<n-.>^(^)  -e( 


*^IN  . ^ITN^ 

"'it  “i  y 


DEN  = 


■ (^)  ■'eN  ^ ^ ^ •'in/”]  • 


ANLP  = 


DIFFE,  = (l/DEN)(kT./m]  ^ - -2!^ 

2 ' I ""it  mj  J 

AMI  p / 1 \ [‘'ITN  ( ’^ItN  / *'lN  ^ITn\' 

“ (m  • DENj  “m][“  \ / ^nij^  mj  ) 

DIFFITj  = (l/DENXkTj/mj^)  + . 

DIFFITj  = (l/DEN)  • 

ANLP  , - '1/DENI  (i?m/m,T)  • 

\ ®/ 


/ ^ Ntt  \ '^Q  *^eN 

DIFFITj  = (l/DEN)  (-5^1 


ANLP 


When  Nj^-  0, 


then  Tj^  -•  0 and 


k(Tg  + Tj)  3Ng 

^e  “*■  ■ _ Tj.  3x 

"'l^IN 


I 

i 


which  is  the  usual  expression  for  amblpolar  diffusion  In  the  presence  of  a single 
positive  Ion  species. 

Substituting  Eqs.  (32)  and  (33)  Into  Eqs.  (26)  and  (27),  and  then  substituting 
Eqs,  (26)  and  (27)  into  Eq.  (25)  gives 


r ^e®‘ 

r— £ = iue  + — 2— 

9x  L “ ■ 


^(lu)  - 7jj^)  N^e^du  - v^j^) 
■ DENOM  Mj  • DfeNOM 


It  Is  reasonable  to  assume  that  the  scale  length  for  electron  density  changes 


^N_-  [ax  J 


Is  long  compared  with  the  laser  wavelength  in  the  plasma  produced  in  front  of 

the  target; 

(55)  » K . 

e ^ 

It  is  convenient  to  rewrite  Eq.  (42)  in  the  form 


■s — = lije  E - (tE  = iue  KE 
8x  ®o  y y O y ’ 


where  <t  the  conductivity,  is  given  by 

a - ‘ ^eN>1 

^ ' DENOM  L m mj[  • 


and  the  dielectric  constant  K is  given  by 


K = 1 


Combining  Eqs.  (42)  and  (24)  yields  the  wave  equation 


3‘^E  2 

-^+— j-KE  =0  . 

dx  c ^ 


Since  the  electron  density  N and  electron  temperature  T_  are  functions  of  the 

“ © 

electric  field  Intensity  Ey,  and  since  K Is  a function  of  N^  and  T^,  Eq.  (44)  is 


nonl’near  in  E . However,  the  use  of  assumption  (55)  permits  Eq.  (44)  to  be 

y 28 

solved  using  the  WKB  perturbation  method  (Ginzburg,  chapter  39): 


E - njm 


(45) 


where  Is  Incident  field  Intensity  at  position  x = x^. 

Thus,  under  these  55  assumptions,  the  transport  equations  coupled  with 
Maxwell's  equations  reduce  to  Eqs.  (7),  (8),  (9),  (10),  (11),  (12),  (16),  (17),  (18), 
(20),  (21),  (22),  (23),  (24),  (34),  (36),  (40),  (41),  and  (45).  It  should  be  noted  that 
the  terms 


can  be  deleted  from  Eq.  (20),  because  it  was  shown  in  reference  9 that 


Vjj*  « t3kTg/m]^/^  . 


As  indicated  In  reference  9,  the  electron  temperature  T^,  the  electron  flux 
density  r^,  and  indirectly,  the  electron  density  N^,  ionization  frequency  Vj,  elec- 
tron thermal  conductivity  K^,  target  temperature  T,j,  (x  = L),  target  vapor  density 
N,p,  and  so  forth,  are  all  functions  of  the  electromagnetic  field  intensity  | E^| . 
This  implies  that  all  particle  densities  and  temperatures  are  mutually  coupled. 
The  system  of  nonlinear,  partial  differential  equations  is  solved  by  converting 

them  to  a system  of  partially  implicit,  partially  explicit  (PI  - PE)  difference 
21 

equations. 

The  boundary  conditions  are  such  that  the  target  temperature  T,p  (x  = L)  is 
given  by  Eq.  (1),  but  T,p  (x  = L)  cannot  exceed  T'  given  by  Eq.  (5).  The  target 
vapor  density  N^,  at  the  target  surface  x = L,  is  given  by  the  expression 


t'=t 


t'  = 0 


dt' 


3^  N„(x  = L) 

Nt3  I'q  exp  (-LM/N^j  kT.^  (x  = U)  + 


8x‘ 


(46) 


where  N,j,^  = target  material  density  at  target  surface.  The  momentum  and  energy 
transfer  terms  in  the  transport  equations  are  set  up  in  such  a fashion  that  the  total 

28.  Ginzburg,  V.  L.  (1964)  The  Propagation  of  Electromagnetic  Waves  in  Plasmas. 
Pergamon  Press. 
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momentum  and  energy  of  the  entire  system  is  conserved.  Also,  the  gain  and  loss 
terms  in  the  continuity  equations  are  such  that  the  total  particle  density  Is  con- 
served. Thus,  there  is  no  need  to  impose  particle,  momentum,  and  energy  con- 
servation conditions  at  the  boundaries  of  the  system.  The  boundary  conditions  on 
the  dependent  variables  N^,  Nj^,  Np  T^,  Tj,  and  Tj^  are  free. 

The  electromagnetic  field  quantities  and  are  specified  at  x = -oo  (laser  wave 
source).  The  solution  to  the  laser-target  interaction  problem  is  determined  with 
the  specification  of  the  initial  conditions  on  die  dependent  variables  N^,  , 

Nt.  Tjj.  T^.  Tj.  Tjj.  and  Tpfpr  -oe<  x<  +»  . 

In  this  report,  the  background  gas  is  taken  to  be  xenon.  The  transport 
jfficients,  reaction  rate  coefficients,  and  cross  sections  have  been  selected 
from  a number  of  references.  A list  of  the  physical  constants  and  transport 
coefficients  follows: 

Vj,  = ionization  frequency  for  electrons  impacting  on  the  neutral  particles  in 
the  ground  state, 

V,  8kTg/(»m)]  exp  (-Ej/kT^)  , 

Ej  - ionization  energy  of  ambient  gas  atoms  (xenon) 

= 19.408  X 10"^®  joules, 

■ Bohr  radius  = 0.  53  X lO”^®  m. 

29 

This  resslon  for  Vj  was  taken  from  Lin  and  Teare. 

= ionization  frequency  due  to  electron  impact  on  target  vapor  atoms, 

= 8kTg/(xm)]  (ira^)  exp  (-Ej.p/kTg)  . 

Ej„  ' ionization  energy  of  target  vapor, 

-1  Q 

= 5.  984  X 1.  6 X 10  joules. 

^NNI  ~ ^o^^zation  frequency  due  to  a neutral  atom  colliding  with  an  ambient 
gas  atom, 

= . exp  (-Ej/kTjj)  . 

This  expression  for  was  taken  from  Zeldovich  and  Raizer. 

v^Nit  ~ ionization  frequency  due  to  a neutral  atom  colliding  with  a target 
vapor  atom, 

= lO"^^  [skT.p/(xm.j.)]^/2  . |-Ej,j,/(kT^)  + ij  • exp  j^-Ej,j./kT.pJ  . 

29.  Lin.  S.  C.  and  Tears.  J.D.  (1962)  Avco  Research  Report  115.  Contract  No. 

AF19(604)-7458.  

30.  Zeldovich,  Ya.  B.  and  Raizer,  Yu.  P.  (1966)  Physics  of  Shock  Waves  and 

High-Temperature  Hydrodynamic  Phenomena,  Vol,  1,  Academic  Bress. 
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Vj®*  = Ionization  frequency  due  to  electron  impact  on  an  excited  state  atom, 

= [8kT^/(xm)l^/^(^aj*)  • exp  - {(Ej  -E^^)/kT^]. 

Egj^  = excited  state  energy, 

= 14. 4 X lO"^®  joules. 

- total  recombination  rale  coefficient  for  radiative  recombination  and 

electron-electron-ion  recombination  (from  analogous  calculations 

31 

performed  for  hydrogen,  as  given  in  McDaniel  ). 

q = photoionization  rate  for  ionization  of  ambient  gas  atoms  by  continuum 
radiation, 

where 

y = hi//(kTj^). 

V = frequency  of  continuum  radiation, 

= spectral  intensity  of  continuum  radiation, 

= cross  section  for  the  absorption  of  a photon  hv  by  an  atom  in  the  nth 
state  and  removal  of  an  electron. 

30 

This  expression  for  q is  given  by  Zeldovich  and  Ralzer.  For  a hydrogen- 
like  atom  whose  remainder  charge  is  Z in  the  nth  quantum  state 

-99  9 

= 7.9  X 10  ^ (n/Z , in  meters  squared  , 

where  is  the  minimum  frequency  of  a photon  still  capable  of  removing  tin  elec- 
tron from  the  nth  level  (Kramer's  formula). 

For  the  steady  state  case,  where  the  temperature,  density  distribution,  and 
radiation  field  are  in  local  thermodynamic  equilibrium,  the  spectral  intensity  is 
given  by 


X 

yx'=/  Kj^I^pexp 

X 

-/  Kj^dx" 

d*'  + lyo  ®*P 

X 

-f  K|,dx" 

’'o 

X' 

*0 

31.  McDaniel,  E.W.  (1964)  Collision  Phenomena  in  Ionized  Gases,  John  Wiley 
and  Sons. 
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where 


VO 

Kl. 


K!. 


K 

i/p 


- Intensity  of  any  Incident  external  radiation, 

= absorption  coefficient  lowered  by  Induced  emission, 
*K^[1  -exp(-hi//kTjj)l  . 

= absorption  coefflclont, 

= Planck  distribution  for  radiative  equilibrium. 


= (2hv^/c*) 


’ 1 

_ exp  (hv/kTj^)“^_ 


For  continuous  absorption  of  light  In  inert  monatomic  gases  when  hi/>  Ep  an 
approximate  expression  for  Is  given  by 


in  Inverse  meters  , 


where 

Vl  = . 

y =hv/kTpf  . 

When  the  plasma  In  front  of  the  target  is  sufficiently  hot  to  radiate  appreciably, 
then  1 (the  plasma  is  optically  thick  to  the  continuum  radiation).  In  this 

case, 

(i^r) ’.i  0 ■ . 


(i) 


7.9  (10"^^) 


ifi  . r (1  - e-y> 

(he)  •{  (e^  + 


(e^  + l)y 


dy  . 


-22  2 Ej  _ ^1 

s 4(l/h)  7.9  (10 

(hc)^  y 


Qex  ~ photo-lonizatlon  rate  for  Ionization  of  excited  state  atoms  by  continuum 


radiation. 


= 8(l/h)  7,9  (10‘^^) 
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While  the  continuum  radiation  Is  such  that  the  plasma  appears  optically  thick, 
the  resonant  line  radiation  can  diffuse  because  of  the  finite  line-width.  The  follow- 
ing expressions  describing  the  transport  processes  of  the  "wings"  of  the  resonant 

25 

radiation  line  are  taken  from  Myshenkov  and  Ralzer. 


Drad  ■ < 

- effective  resonant  radiation  diffusion  coefficient, 

i E mean-free  path  of  resonant  radiation  quanta, 

- 0.7901  X 10"'^  m, 

T = average  lifetime  of  excited  state, 

=■  3.74  X 10‘®  sec. 

T = average  time  for  resonant  radiation  to  escape  from  plasma  (laterally, 
In  y - z plane), 

R = average  distance  from  center  of  plasma  to  lateral  boundary  (approx. 
0.01  m). 

= collision  frequency  for  production  of  excited  state  atoms  due  to 
electron  Impact  on  ambient  gas  atoms. 


■/ 


“.x'o*' 


. 10-22(2/^  • ««P  (-S.  . l/s,|  . 

f^  = Maxwellian  electron  distribution  function, 

= exp  (-mv^/(2kT^))  . 

'^ex  = V 2E^x"/-^  • 

V = electron  velocity. 

Ogjj  = cross  section  for  excitation, 

= 10“^^  (from  Massey  et  al, 

= m/(2kTg). 

32.  Massey,  H.  S.W. , Burhop,  E.  H.  S. , and  Gllbody,  H.B.  (1969)  Electronic 
and  Ionic  Impact  Phenomena.  Oxford  Clarendon  Press. 


* electron-neutral  collialon  frequency  for  momentum  transfer^ 

* Nj^4ir(m^/airkT^)^^*  f v®  exp  {-(J^  v^l  dv, 

0 

^eN  ' momentum  trenafer  croas  section, 

* 3(10"^®)  m®  (from  Massey  et  al.®®). 

Vej  ■ electron-ion  collision  frequency, 

N (lnA)T  - .. 

» - OTST (10  ®)  (from  Holt  and  Haakeir®). 

A =<d/Pc- 

= 5.55(10“®) 

ijj  ^ 90.C  (T^/N^)^/®. 

*'lN  ~ ion-neutral  collialon  frequency, 

’ **N®IN'^r 

vj  . |3kTj/mjl 

*^IN  * cross  section  for  momentum  transfer, 

* e/  (m jp  Njj  Vj) . 

Mjjj  « mobility, 

« 0.  6(10’®)  m®/(volt-sec)  (from  McOaniel®^. 

- collisicn  frequency  for  electron-excited  state  atoms, 

erx 

'^lex  ~ ion-excited  state  collialon  frequency, 

W^ex/^N^’ 

'^ITN  ~ i*rget  vapor  ion-neutral  gas  collision  frequency. 

’ Nn^inV__ 

VjT  ' VsJsTj/mj^  . 

*'nTN  ‘ target  vapor— neutral  gas  collision  frequency, 

Vjjip  = ^ 3kTij,/m,j,  . 
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'S' 


« thermal  conductivity  of  ambUnt  |aa, 

> thermal  conductivity  of  target  vapor, 
Ml/8)(8kT^vm,j.)'/^  . 

^(4,aJ> 


■ electron  thermal  conductivity. 


K, 


m„ 


'eN 


'el 


'IN 


5N  k^T 

~mi>  *"  Shkarofaky  et  al. 


el 

< ion  thermal  conductivity, 

SN.k^T 
Jf L. 


(from  Shkarofaky  et  al. 


- target  atom  mesa, 

= 4.  514  X 10“^®  kg  (aluminum). 

< apecific  heat  of  target  material, 

» 0.  9 X 10®  joules/(kg  - ®K)  (aluminum). 

= heat  of  vaporization, 

» 1.076  X lo''  joulea/kg. 

> density  of  aluminum, 

= 5.986  X 10^®  kg/m®, 

= relative  fractional  energy  loss  for  elastic  electron  — neutral  collisions, 

= 2m^/mj^. 

= G. 


'eN  • 

(mj  + mj^)^ 


■=  1/2. 


« rate  of  emission  of  target  vapor  at  target  surface. 
= N.p(x  * L)  exp  l-LM/N^kT.j.(x  = L)J. 


33.  Shkarofaky,  I,  P.,  Johnston,  T.W. , and  Bachynskl,  M.  P.  (1966)  The  Particle 
Kinetics  of  Plasmas,  Addison-Wesley. 
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M - *tomte,w*lfht  of  Urgot  mMorUl. 

■ Avogkdro'i  nutnbor. 

Shock  formation  In  the  ambiant  gas  la  doacribad  by  adding  a paaudo-vlaeoua 
praaaura  tarm  (QSV)  to  Eqa.  (11),  (17),  and  (23)  (Rlchtmyar  and  Morton^S: 

In  Eq.  (17), 

In  Eq.  (23), 

- <rN/ ^ (^)  <^'n/ V * 

whara 

QSV  . 2(Ax)*(mj^/Njj) 

Hera,  Ax  « apatlal  atap  slaa.  The  addition  of  tha  paaudo-vlaeoua  praaaura  QSV 
allowa  N^,  and  to  vary  amoothly  ovar  a ahoekad  raglon  with  a ahock-front 
thlcknaaa  about  aqual  to  aavaral  Ax. 

a 

It  waa  axplalnad  In  a pravloua  raport  that  tha  numarlcal  Intagratlon  of 
Eqa.  (7),  («).  («),  (10),  (ll),  (12),  (16),  (17),  (18),  (20).  (21),  (22).  (23),  (24), 
(34),  (36),  (40),  (41),  and  (45)  la  graatly  facllltatad  by  braaklng  tham  down  Into 
aubgroupa.  Flrat,  tha  partial  dlffarantlal  aquatlona  ara  convartad  to  a sat  of 
partially  Implicit,  partially  explicit  (PI  - PE)  dlffaranca  aquatlons.^^  If  thaaa 
aquations  ara  proparly  grouped,  than  tha  PI  • PE  dlffaranca  aquations  within  each 
group  may  be  solved  Indapandantly  of  tha  aquatlona  In  other  groupa. 

Tha  system  of  PI  - PE  difference  aquatlona  may  be  written  In  matrix  form  an 

Aa-B  . (47) 

where  £ Is  a column  va'^tor  representing  tha  unknown  dependent  variables  at  the 
new  time  atap  t At,  * matrix  composed  of  known  transport  eoafflclanta  at 
the  old  time  t,  and  £ la  a column  vector  composed  of  known  dependent  variables 
and  transport  coefficients  at  tha  old  time  t.  Tha  determination  of  tha  dependant 
variables  (Njj,  N^,  Np  N,j.,  Nj^,  Fj^,  r^,  etcetera)  at  tha  new  time  atap 
t At  Involves  tha  inversion  of  tha  matrix  Tha  determination  ot  la  graatly 
facilitated  by  grouping  the  aquations  as  follows: 

Group  I,  ambient  gas  atoms  Eqs.  (11).  (17).  and  (23); 

Group  n,  excited  state  atoms  -•  Ec.j.  (8),  and  (18); 
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Group  FI,  electrons  and  iona ->  Eqs.  (7),  <10),  (20),  (21),  (34),  (36),  (40), 
and  (41); 

Group  IV,  target  vapor  Eqs.  (12).  (16),  and  (22). 

By  following  this  particular  subdivision  into  groups,  the  matrix  A to  be 

ss 

inverted  for  each  subgroup  is  much  smaller  than  the  matrix  corresponding  to  all 
equations  considered  simultaneously.  This  subdivision  into  groups  results  in 
great  saving  of  computer  storage  requirements  and  computer  running  time. 

A subroutine  called  Define  has  been  written  to  create  the  appropriate  entries 
in  and  B for  each  subgroup  of  equations.  Since  the  original  differential  equations 
were  of  second  order  in  the  spatial  variable  x,  the  corresponding  difference  equa- 
tions are  of  second  order  in  the  spatial  step  size  Ax.  This  implies  that  A is  a band 

fS 

structured  matrix,  with  many  co-diegmals  identically  zero.  The  band  structure 

of  A is  decomposed  in  Define  into  a sic.glr  column  vector,  as  described  in  the  IBM 

® 34 

scientific  subroutine  package  entitled  GELB  (Gauss  Elimination  Band).  This 

scientific  subroutine  package  GELB  solves  M simultaneous  equations  in  M unknowns 

35 

very  efficiently  by  Gauss  elimination  with  column  pivoting  only  (Hildebrand  ). 

The  high  energy  laser-target  Interactions  computer  program  consists  of  a 
main  program  that  calls  about  40  subroutines.  The  initial  conditions  for  the 
particle  densities,  fluxes,  and  temperatures  are  read  in  as  data.  The  value  of  the 
initial  target  temperature  is  specified,  and  the  value  of  the  incident  laser  field 
intensity  is  also  specified.  Parameters  determining  the  degree  of  implicit- 

ness of  the  diiference  equations  and  the  spatial  step  size  Ax  are  initially  given. 

The  electric  field  E and  magnetic  field  H are  determined  at  each  time  step 
y z 

from  Eqs.  (45)  and  (24). 

The  step-by-step  time  integration  of  these  transport  equations  is  accomplished 
by  calling  the  subroutines  Define  and  Gelb  from  the  main  program  for  each  sub- 
group of  equations.  This  determines  all  the  dependent  variables  along  the  new 
time  line  (t  + At).  The  time  step  At  is  automatically  adjusted  in  the  main  program 
by  requiring  that  the  successive  changes  in  the  dependent  variables  lie  between 
certain  preselected  tolerances.  Let 

U^(t  + At)  - U,  (t) 

"^°^k  ^ Uj^(t  + At)  ■ 

where  Uj^  is  any  dependent  variable  (such  as  N^,  and  so  forth).  If  TOLj^  for 
all  variables  is  less  than.TOLMIN,  the  time  step  At  is  increased  by  ?0  percent. 


34.  IBM  Application  Program  (1970)  System/360  Scientific  Subroutine  Package, 

Version  m,  edition  GH20-0205-4.  

35.  Hildebrand,  F.  B.  (1056)  Introduction  to  Numerical  Analysis,  McGraw-Hill. 
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To  reduce  the  truncation  error,  the  time  atep  At  la  cut  In  half  If  TOLj^  la  greater 
than  TOLMAX  for  any  dependent  variable.  Alao,  to  Increaae  the  efficiency  of  the 
program.  If  TOL^  la  lea  a than  0. 1 TOLMIN  for  a particular  aubgroup  of  equatlona 
for  ten  conaecutlve  time  atepa,  then  that  aubgroup  of  equatlona  la  "frozen  out"  or 
not  Integrated,  for  ninety  conaecutlve  time  atepa.  All  the  transport  coefficients 
such  as  Vj,  Vj  , Vj^,  a^,  and  so  forth,  are  evaluated  at  each  time  step 

In  separate  computer  subroutines. 

The  computer  output  consists  of  listings  of  all  the  dependent  variables  as  a 
function  of  position  at  each  time  step.  Also,  the  laser  field  Intensity  at  the  target 
surface  and  the  target  surface  temperature  are  given  at  each  time  step. 


3.  CONCLUSIONS 

It  was  explained  In  reference  9 how  the  high  energy  laser-target  Interactions 
program  was  checked  out.  Let 

P = degree  of  Implicitness  for  first  order  spatial  derivatives, 

S = degree  of  Implicitness  for  second  order  spatial  derivatives, 

Z = relative  fraction  taken  for  forward  differences  along  the  new  time  line, 

V = relative  fraction  taken  for  forward  differences  along  the  old  time  line. 
For  example,  S = 0 corresponds  to  a fully  explicit  Integration  scheme  for  second 
order  spatial  derivatives  and  S = 1 corresponds  to  a fully  Implicit  scheme.  Also, 
Z = 1 corresponds  to  caking  forward  differences  along  the  new  time  line,  and 
Z = 0 corresponds  to  taking  backward  differences  along  the  new  time  line. 

In  this  report,  the  spatial  step  size  Ax  Is  nonuniform.  The  step  size  Ax  was 
taken  equal  to  2.  676  X 10  meters  at  the  target  surface.  Each  succeeding  step 
size  Is  Increased  3.  5 percent  as  one  recedes  away  from  the  target  surface.  The 
maximum  number  of  Incremental  steps  equals  154.  The  total  distance  considered 
equals  1.  521  mm.  Maximum  numerical  stability  was  achieved  with  the  following 
scheme: 

S = 0.  For  Group  I (neutral  atoms), 

> Group  n (excited  state  atoms). 

P = 0.  5 ) Group  IV  (target  vapor). 

S = l.O-j 

J For  Group  III  (electrons  and  Ions). 

P=  1.0  / 

Let  J represent  an  Index  to  count  the  position  of  the  atep  size  Ax,  so  that 
J = 1 corresponds  to  the  target  surface  and  J = 155  corresponds  to  the  position  of 
the  Incident  laser  beam.  The  parameters  V and  Z are  chosen  so  that 
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V = 0. 5 , 
Z‘0.5  . 


This  corresponds  to  taking  central  dil'ferences. 

In  this  report,  the  laser  frequency  u is  set  equal  to  1.77826  x 10^^  rad/sec 
(CO2  laser  frequency).  The  target  is  aluminum  and  the  ambient  gas  is  xenon.  In 
the  figures,  the  initial  conditions  on  the  variables  Tjj,  Nj^,  T^j,  Tj,  T^, 

Ne,  Njip,  Nj,  and  T,p  are  obtained  by  first  selecting  reasonable  profiles  for  these 
variables  and  then  letting  the  program  run  with  a small  value  for  (incident 

field  intensity)  until  the  successive  changes  in  the  variables  become  relatively 
small.  Each  case  is  illustrated  by  a set  of  two  figures;  one  figure  represents  the 
particle  densities  as  a function  of  position  at  a given  instant  in  time  and  another 
figure  represents  the  particle  temperatures  as  a function  of  position  at  a given 
instant  in  time.  In  these  figures, 

O — > electrons, 

□ — neutral  gas  atoms, 

0 — > target  vapor  atoms, 

^ — > in  the  density  plots,  excited  state  atoms, 

A — ^ in  the  temperature  plots,  ion  temperature. 


The  distance  x is  normalized  in  these  plots:  x/L,  where  L = total  laser  path 
length  - 1.  521  mm. 

In  Figures  1 and  2,  is  set  equal  to  1 volt/meter  (2.  7 X 10*^  watts/cm^) 

and  it  may  be  noted  that,  after  an  elapsed  time  of  0.  257  psec,  all  density  and 
temperature  profiles  are  fairly  uniform  in  space.  The  final  density  and  tempera- 
ture profiles  in  Figures  1 and  2 are  used  as  initial  conditions  in  all  subsequent 
calculations.  In  this  report,  the  incident  field  starts  from  a very  low  value 

(1  volt/meter)  and  rises  as  the  square  root  of  the  elapsed  time.  This  implies  that 
the  incident  flux  rises  linearly  with  time,  and  the  rise  time  is  adjusted  so  that  the 
incident  flux  reaches  its  peak  value  after  an  elapsed  time  of  0. 1 ^sec. 

In  Figures  3 and  4,  Ejj^^  = 6 X 10®  volts/meter  (0.95  X 10^  watts/cm^),  and 
the  elapsed  time  = 0.  618  psec.  Ionization  due  to  a neutral  atom  (ambient  gas  or 
target  vapor)  colliding  with  a target  vapor  atom  is  neglected  (•'jjjgj.j.  = 0).  The 
heat  diffusion  in  the  target  vapor  away  hrom  the  target  surface  may  be  noted  in 
Figiure  4.  In  Figures  5 and  6,  Ejj^^  = 6 X 10®  V/m,  and  the  elapsed  time  = 1. 41 
psec.  Also,  = 0.  Comparing  Figure  6 with  Figure  4,  it  may  be  noted  that 

the  diffusion  of  heat  away  from  the  target  surface  continues  in  the  target  vapor. 
However,  it  may  be  noted  that  there  is  no  initiation  of  a laser  absorption  wave 
(LAW)  at  this  incident  power  density  (0.95  x lo"^  watts/cm^). 
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In  Figures  7 and  8,  * 2X10*^  V/m  (1.06X10^  watts/cm^)  and  the  elapsed 

time  ■ 0.048MS6C.  The  ionization  due  to  a neutral  atom  colliding  with  a target  vapor 
atom  is  neglected  0)*  Figures  9 and  10,  the  conditions  are  the  same  as 

in  Figures  7 and  8,  except  that  the  elapsed  time  >>  0.  103  psec.  It  may  be  observed 
that  the  electron  density  is  increasing  at  a distance  of  about  one  wavelength  from  the 
target  surface.  This  indicates  the  ignition  of  a LAW.  In  Figure  10,  it  may  be  ob- 
served that  the  ambient  gas  temperature  reaches  a peakvalue  of  about  10^  °K,  near 
the  target  surface.  Also,  the  target  vspor  temperature  reaches  a peak  value  of 
about  10  K near  the  target  surface.  This  is  physically  unrealistic,  because  at  this 
temperature  two  target  vapor  atoms  would  collide  and  ionize.  This  shows  that  target 
vapor  atom-target  vapor  atom  impact  ionization  may  not  be  neglected  in  the  ignition 
of  LAW'S  ^ 0)*  target  vapor  reaches  a high  temperature  because  the 

target  vapor  atoms  are  being  rapidly  emitted  from  the  surface  and  are  diffusing  into 
approximately  one  atmosphere  of  xenon  gas,  so  tha*  the  target  vapor  is  undergoing 
rapid,  nearly  adiabatic  compression.  The  target  >or  hasn't  been  able  to  transfer 
all  of  its  excess  energy  to  the  xenon  gas  on  this  time  scale  (0.  103  lusec). 

In  Figures  11  through  18,  = 6 X 10®  V/m  (0.  95  X 10^  watts/cm^).  Also, 

target  vapor  atom  — target  vapor  atom  impact  ionization  is  not  neglected 
^"^NNIT  Figures  3 through  6 with  Figures  11  through  18 

^^INC  “ ® ^ Vim  in  all  these  figures),  it  may  be  noted  that  the  densities  and 
temperatures  are  about  the  same  when  = 0 and  ^ 0 for  comparable 

elapsed  times.  Also,  for  the  case  when  * 0,  it  may  be  noted  that  there  is 

no  ignition  of  LAW'S  at  = 6 x 10®  V/m. 

In  Figures  19  through  24,  Ejjjq  = 2 x 10^  V/m  (1.06  X 10®  watts/cm^).  Tar- 
get vapor  atom —target  vapor  atom  impact  ionization  is  not  neglected  *0). 

The  elapsed  time  in  Figure  23  (0. 074  p sec)  is  roughly  comparable  with  the  elapsed 
time  in  Figure  9 (0. 103  psec).  However,  the  peak  target  vapor  density  in  Figure 
23  is  about  one  order  of  magnitude  less  than  in  Figure  9.  Also,  it  may  be  seen 
that  the  peak  target  vapor  temperature  is  almost  an  order  of  magnitude  less  in 
Figure  24  compared  with  Figure  10.  One  may  discern  from  Figure  23  the  ignition 
of  a LAW  by  an  increase  in  the  electron  density  profile  at  a distance  of  about  two 
laser  wavelengths  from  the  aluminum  target  surface. 

In  Figures  25  through  28,  = 2 x 10^  V/m  and  ^0,  just  as  in 

Figures  19  throu^  24,  except  that  in  Figures  25  through  28,  the  diffusion  of  the 
electrons  and  ambient  gas  ions  is  uncoupled  from  the  diffusion  of  the  target  vapor 
ions.  This  uncoupling  is  achieved  by  letting  DIFFE2  — 0 in  Eqs.  (40),  DIFFI2  — 0 
and  ANLPI  — 0 in  Eq.  (41),  and  (Nj,j,/N^)  — 0 in  the  expressions  for  DIFFEl  and 
DEN  under  Eq.  (41).  Comparing  Figures  25  and  26  with  Figures  19  and  20,  and 
Figures  27  and  28  with  Figures  23  and  24,  it  may  be  observed  that  this  uncoupling 
procedure  produces  no  discernible  effects  in  the  density  or  temperature  profiles 
for  comparable  elapsed  times. 
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In  Figures  29  through  38,  Ejj^^  * 2 X 10*^  V/m  (1.06  X 10®  watts /cm^)  and 
*'nnIT  * ®*  marks  In  the  density  and  temperature  profiles  at  x > 0. 9L 

Indicate  a scale  change  in  Figures  29  through  38.  In  Figures  29  and  30,  the 
elapsed  time  = 0. 079  fisec.  The  Ignition  of  a LAW  may  be  clearly  observed  in 

Figure  29  by  the  Increase  in  the  electron  density  profile  at  a distance  of  about 

-5 

10  meters  (one  laser  wavelength)  from  the  target  surface.  In  Figure  29,  it  may 
be  seen  that  the  aluminum  vapor  layer  extends  from  x = 0.  94L  up  to  x = L. 

In  Figure  30,  it  may  be  noted  that  the  target  vapor  temperature  reaches  a 
value  of  about  2 X 10®°K  at  the  outer  edge  of  the  aluminum  vapor  layer  (x  = 0.  94L), 
due  to  nearly  adiabatic  compression.  The  vapor  emerging  at  the  target  surface  la 
about  2000°K. 

From  Figure  29,  it  is  seen  that  the  target  vapor  density  N,p  increases  as  x 
increases  towards  the  target  surface  (x  = L).  Also,  from  Figure  30  it  is  seen  that 
the  target  vapor  temperature  T,p  increases  rapidly  as  the  edge  of  the  vapor  layer 
is  reached,  and  then  decreases  more  slowly  as  x approaches  the  target  surface. 
Since  the  rate  of  production  of  electrons  by  the  process  A1  + A1  -•  A1  + Al^  + e 
increases  as  N,p  and  T,p  increase,  it  is  not  surprising  that  the  electron  density 
buildup  shown  in  Figure  29  occurs  somewhere  in  the  middle  of  the  vapor  layer 
extending  from  x = 0. 94L  to  x = L.  Figures  31  through  38  show  the  farther  develop 
ment  of  the  laser  absorption  wave.  In  Figure  35  (elapsed  time  = 0. 1 psec)  and  in 
Figure  37  (elapsed  time  = 0. 102  psec),  the  target  vapor  density  exceeds  the  xenon 
gas  density.  Also,  from  Figures  36  and  38  it  may  be  seen  that  close  to  the  target 
surface  the  xenon  gas  temperature  T^  has  increased  to  over  2000*^K,  having  been 
heated  primarily  through  elastic  collisions  with  the  aluminum  vapor  atoms.  Now, 
the  xenon  gas  is  sufficiently  hot  so  that  the  process  Xe  + Xe  Xe  + Xe^  + e also 
plays  a role  in  the  electron  density  buildup.  This  process  accounts  for  the  "spike" 
on  the  electron  density  profile  near  the  target  surface  that  occurs  In  Figures  35 
and  37. 

To  summarize,  the  early  stages  in  the  ignition  of  a xJlW  from  an  aluminum 

8 2 

target  in  xenon  were  observed  at  a power  density  of  10  watts /cm  . Starting  with 
a low  density  of  "priming"  electrons  of  l0®/cm®  (due  to  "prompt"  field  emission 
in  time  scales  less  than  or  of  the  order  of  1 nsec),  the  first  mechanism  responsi- 
ble for  substantial  electron  density  buildup  is  A1  + A1  — A1  + Ai  + e.  Later,  as 
the  ambient  xenon  gas  is  heated  by  the  hot  target  vapor  through  elastic  collisions, 
a second  electron  density  buildup  mechanism  starts  to  play  a role: 

Xe  + Xe  — Xe  + Xe'*'  + e. 

Shock  waves  can  be  generated  in  the  ambient  gas,  but  not  in  the  electron  gas, 
since  v ® is  always  much  less  than  ^3kT_/m  . Also,  there  is  considerable 

X O' 
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difference  between  the  temperatures  of  the  electrons,  ions,  target  vapor  atoms, 
and  xenon  gas  atoms.  The  computer  program  that  has  been  developed  can  be  used 
to  study  the  further  development  of  LAW'S.  Also  the  Ignition  and  propagation  of 
LAW'S  from  various  materials  (such  as  dielectrics,  or  coated  metal  targets)  may 
be  investigated. 


45 


^ 10>9 

in 

Z n 
UJ  10*’ 
Q 


10" 

0. 


Figure  1.  Particle  Density  vs  Normalized  Distance 
(F  « 0 MV/m;  Time  = 0. 257  X 10"6  sec) 
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NORMALIZED  DISTANCE 

Figure  2.  Temperature  (°K)  vs  Normalized  Distance 
(E  = 0 MV/m;  Time  = 0. 257  x 10-»  sec) 
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ZMPERfiTURE 
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Figure  7.  Particle  Density  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.048  X lO"®  sec) 
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Figure  8.  Temperature  (°K)  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.048  X 10"“  sec) 
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't  Figure  20.  Temperature  (°K)  vs  Normalized  Distance 

' I (E  = 20  MV/ms  Time  = 0. 046  x 10"®  sec) 
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Figure  21.  Particle  Density  vs  Normalized  Distance 
(E  = 20  MV/mj  Time  = 0.061  X lO'S  sec) 
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Figure  22.  Temperature  (^K)  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.061  X 10-8  gee) 
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Figure  23.  Particle  Density  vs  Normalized  Distance 
(D  = 20  MV/m;  Time  = 0.074  X 10~6  sec) 


Figure  24.  Temperature  (°K)  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0. 074  X 10“®  sec) 
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Figure  25.  Particle  Density  vs  Normalized  Distance 
(B  = 20  MV/m;  Time  0.046  x 10’“  sec;  uncoupled 
diffusion) 


Figure  26.  Temperature  (°K)  vs  Normalized  Distance 
(E  » 20  MV/m;  Time  » 0.046  x 10’®  sec;  uncoupled 
diffusion) 
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igure  27.  Particle  Density  vs  Normalized  Distance 
2 » 20  MV/m;  Time  = 0.074  X 10"®  sec;  uncoupled 
Lf fusion) 
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Figure  28.  Temperature  (°K)  vs  Normalized  Distance 
(E  » 20  MV/m;  Time  = 0.074  X 10"®  sec;  uncoupled 
diffusion) 
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Figure  29.  Particle  Density  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.079  X 10’®  sec) 


0 0.90  0.95  1.00 

NORMALIZED  DISTANCE 

Figure  30.  Temperature  (°K)  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  » 0. 070  X 10"®  sec) 
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Figure  31.  Particle  Density  vs  Normalised  Distance 
(E  = 20  MV/mj  Time  = 0. 080  X lO*®  sec) 
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Figure  32,  Temperature  (°K)  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.080  X 10’®  sec) 
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Figure  33.  Particle  Density  vs  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.093  X i0“®  sec) 
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Figure  34.  Temperature  (°K)  va  Normalized  Distance 
(E  = 20  MV/m;  Time  = 0.083  X lO'®  sec) 
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